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It is shown that fundamental solutions K'^{x,y;t) = {a;|e“'^ *|y) of the non-stationary 
Schrodinger equation (Green functions, or propagators) for the rational extensions of the Harmonic 
oscillator = Hose + Al/'^ are expressed in terms of elementary functions only. An algorithm to 
calculate explicitly Ka for an arbitrary cr G N is given, compact expressions for and are 

presented. A generalization of the Mehler’s formula to the case of exceptional Hermite polynomials 
is given. 


Introduction 


In this work we present new examples of exactly-solved propagators in one-dimensional quantnm mechanics. We 
consider rationally extended Harmonic oscillators [l|. In this case the evolution of wave packets is periodic, we have 
so called isochronous anharmonic oscillators 0. This periodicity is related with the quasi-equidistant structure of the 
spectrum of the Hamiltonian Q . 

The simplest rational extension is given by the potential Q 
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which leads to the quasi-equidistant spectrum Q, = n + where n G No \ {1, 2}, for the stationary Schrodinger 
equation. 

Another example is the two-well perturbation of the oscillator 


y{2.3}[^] = ^+2 
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with the quasi-equidistant spectrum, i?„ = n + i, n G No\{2,3}. 

Note that each rational extension of the Harmionic oscillator is defined by a sequence of levels tr G N which are 
deleted from the spectrum by Darboux-Crum transformations [I| . Darboux transformations represent a powerful tool 
to manipulate physical properties of one-dimensional quantum systems, In the case of shape-invariant potentials, 
relations between propagators were studied in [ 1 , 0 ] • A more general approach for the calculations of propagators for 
potentials generated by an arbitrary chain of Darboux transformations were developed in (MU- 

As a particular example, propagators for the family were defined through a generating function 

S{x,y,t\J) which contains the error-function m- Here we will further simplify and extend this result. First, we will 
show that the propagator K'^ for an arbitrary potential is expressed by elementary functions only. In the case of 
potentials © and (HI) we get the following propagators 


A:^^’2>(x, y;t) =e "^'^Koseix, 


Ar^2.3}(2., y-t)=e '^^^Koseix, y;t) il- 


4i sin t [xy — e' *] \ 

(l+x2)(l+y2) j ’ 

8i sint [xy{x^y'^ ~ 3) — 3(^2 + y'^) cost — 3i {x'^y'^ + l)sint] 
(3 + x4)(3 + y4) 


where the propagator of the Harmonic oscillator [1^ is used 
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Second, we will define a rational anzatz to compute propagators K'^. In the general case, propagators for the 
rationally extended oscillators have the following structure 


^[[-i]]+i 

E Q^(x,2/)e-‘'=* 


= Kosc{x,y-t)- 
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E Qlix,y) 
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where Q'j^{x,y) = Q'^{y,x) are some polynomials that can be determined iteratively, which is more efficient than 
the method based on the generating function [HI- The polynomials QZ(x, y) allow also calculate Green functions 
Ga-{x,y, E) and generalize Melher’s formula for the x-Hermite polynomials |1^ . 

I. HARMONIC OSCILLATOR, HERMITE POLYNOMIALS AND POTENTIALS WITH 

QUASI-EQUIDISTANT SPECTRUM 


Consider the hamiltonian of the Harmonic oscillator 

r.2 
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■^xx 
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with eigen-functions defined in terms of probabilistic Hermite polynomials 
^|;n{x) = Pnilen{x)e~^ , P„=(n!v27rj 

Rational extensions are defined as the following perturbations of the Harmonic oscillator [If 

where a, 

cr = {fci,fci + 1,... + 1} , |cr|=2M, 


(7) 


( 8 ) 


(9) 
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is a strictly increasing sequence of natural numbers, or a Krein-Adle sequence (El |. 

Following standard notations of the Mathematica program language we denote by cr[[i]] and by (t[[— 1]] f-th element 
and the last element of this sequence, respectively. A sequence of natural numbers appearing as a subscript, for 
instance, ij^a-ix), implies a set of elements, 

■Ip^ix) = = {'^l^nAx),1pn2{x), ■ ■ ■ ,'llJn2M{x)} , 

that is, if {An)n>o is a sequence of elements, then A^ni,...,nm} is a set of elements. This agreement allows us to write 
Wronskians in a compact form 


Wr[i{ja{x),x] = 


i^aiii]]{x) 'ipyi 2 ]]{x) ... i/'a[[-i]](a;) 

V'a[[i]](2;) ■■■ V'(,[[_i]](a;) 


( 11 ) 


Hamiltonians iLosc and H'^ can be embedded into a polinomial SUSY algebra I, El 

2M 2M 

LHoec = H'^L , L+L = n(17osc - a[[j]]), LL+ = ’ 


( 12 ) 


where L is a differential operator of 2M-th order El 

,, , ^ WT[jp^{x)U{f{x)},x] 

Wr[tjja{x),x] 

Operator L maps the oscillator eigen-functions ([5]) to the eigen-functions of the rationally extended oscillator 
t/j^ix) = NnL^n{x) , 


(13) 


( 14 ) 
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where a normalization factor is taken into account 

_ 1 

2M \ ^ 

n (n-cr[[j]]) 1 

0 , n ^ a . 

Using explicit form of oscillator eigen-functions (jS]) and the following identities 

2M 

M T—r 

WT:[%l^a{x),x\ =e 2-Wr[He<,(a;),a;] J_J_Pa[[n]] 

n—1 

^ 2.2 2M 

lnWr[V'CT(a:),a:] =- - -^ In Pcr[[n]] + In Wr[He^(a;), a;] 

n—1 

we can express rationally extended Harmonic oscillators through the Wronskian of probabilistic Hermite polynomials 
only 

Ha = -dl^ + Y “ 2a2^(lnWr[He„(a;),a;]) -b2M. (16) 

Note (see for instance i) that Wr[Hecr(x), x] (for the chosen class of a ) is the polynomial of and 

2M 

deg^Wr[He^(a:),x] = y](cr[[n]] -n + l). 

n—1 


It is convenient to introduce normalized polynomials hn{x) 


hn(x) — p7T,He7T,{x), 


(17) 

and the corresponding normalized exceptional Hermite polynomials 



/i^(x) = iV„Wr[ft,o. U /i„, x] 


(18) 

Finally, we define compact notations 



lU(x) = y^v[ipa{x),x\ , W{x) = Wr[HeCT(x),x] , 


(19) 

Wn{x) = Wr[V',T\{<T[[n]]}(a:),a:], Wn{x) = Wr[He<^\{<^[[„]]}(x),x], 


(20) 

^ Wr[He^(x) U {/},x] 

Wr[Hecr(x), x] 


(21) 

Note that 



2M 
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M '' T —W (2Af — l)a; ^ ^ 

TU(x) = e 2 W{x) [[Pa[[n]], Wn{x)=e ^ 

If P<^[[i]] ’ 

(22) 

n—1 

1=1 


1 

II 

a> 
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II. PROPAGATORS OF RATIONALLY EXTENDED HARMONIC OSCILLATORS 

A. Generating function formalism 



The Schrodinger equation for the Green function reads 
{idt - H)^K(x,y;t) =0, K{x,y,0) = S{x - y). 


( 24 ) 
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If two Hamiltonians Hq and are related by 7V-th order Darboux transformation which remove N levels from 
the spectrum of Hq, then corresponding propagators Kq and Kj^ are related as follows [m, 


KN{x,y;t) = J Kq{x, z]t)u„{z)dz . 


(25) 


In the case of rationally extended Harmonic oscillators b = oo, N = |cr|, and the transformation solutions coincide 
with Harmonic oscillator eigenfunctions Un = 

Using (1221) we can replace Wronskians of wave functions by Wronskians of Hermite polynomials 


2M 


K'^{x,y]t) = ^(-I) 




e — —L, 

OO 


W{y) Jy 


Kosc{x,z;t)e Hecr[[„]](-2:)dz 


(26) 


The occurring integrals Kosc{x, z,t)e Y{en{z)dz can be represented as derivatives of the generating function 
with respect to the auxiliary current J 


j osc {x ■; ^ Hecr[[n]] (-^)^-^ — J=0 — 

dy 

where hm,k are coefficients of the He^- The generating function reads 


a[[n\\ 

k—Q 


cr[[n]],fc- 


d^S{J) 

dJ^ 


,7=0 


1 

S{J\x,y,t) = 2 ® ^ 


sin te 2 


1 xe 2 


i •\/2i sint 


E[J, X, y, t] 


where 


R 


i Jv^ sin te 2 


t xe 2 


i •\/2i sint 


= exp (J(i J sint + x)exp {—it )), 


E[J, X, y, t] = 1 + erf 


T r -—:— ^ iv^ / ii 

— Jv 1 smte 2 -; - ye 2 — xe 2 

2\/shU V / 


Z 

erf (z) = f exp {—t^)dt. 

y/'x J 


Function R is the generating function of rescaled Hermite polynomials (see also Appendix A) 

-2 ' 


R[zy/a, = exp [ xz - ^ 
Va V 2 


n—0 


. z 
n\ 


■‘ 2 * v^Tshit = ^ (1 — e . 
V 2 


^ =^HeH(x): 

For the compact writing we define 
a = —2i sin te' * = I — e , 

In what follows we need derivatives of generating functions R, E and 

Eo[x,y,t] = E[0,x,y,t] = ^1 + erf 
1. J-derivatives of R-function 


iv^ / _ii i*\ \ 

— , xe 2 _ -ye 2 . 

2yihriv ^ J\J 




{k — m)\ 


k\ 


(27) 


(28) 


rHeiri„ (x)e-‘= e-‘HeL“] (x) 


Va'J J^o {k-m)\ ^ dx™ 


(29) 
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2. J-derivatives of E-function 


V 5 J -+1 


Kosc(x,y\t)e e 2 * 2 isint'^ 


i-0 


(x, 2 /)e 




3. x-derivatives of Eq- function 


(30) 


^ (31) 

1=0 

4. Mixed derivatives of E-funetion 


d>^+^Eo[x,y,t] 


= Kosc{x,y;t)e' 


-it 

4 e 2 


1 


2i sint 


i9a;^ \ 9J™+i Jj^q Isisinl 


k-1 


m-\-k 


e Koscix,y;t)e 4 ’' ^ gj,m+fc+i(a::, 2 /)e . 

1=0 


(32) 


In the above expressions qij and Wk^m are some polynomials. 
Now we will substitute these derivatives to calculate 


[He„[[j]] (dj) 5'( J)] = 


-it 3 ;^ 

0 2 4 


^[b1] ^ / \ 

E E Cr 9^-i?[jySe-'‘,^] (a7£;[J,x,y,t]) 

fe=0 m=0 ^ / 


• (33) 


,7=0 


Consider first the following double sum 


•yim 


E ^'^[[1]].*: E 


k\ 


fe =0 


—' 77l!(fc — to)! 
m—Q ^ 


(a7£;[J,x,y,t]) 


,7=0 


We will change the order of summation (k.m —1 m.k), that is we will fix to and calculate first sum by fc G 
(cr[[j]] - TO,cr[[j]]). 


^ 1 nm 

E E (X) 07 E[J,x.i,,il) 

m=0 


J J =0 


<^IW . <^[[1]] ™ 

E = E =rH4"'(-) 


m—0 


k—m 


J J =0 


In the last expression we can change the inferior limit k = m to k = 0 since (x) = 0 when k < m, 

<J]] 


, -i<^[[l]]* 


'^[[11] 

E “1 ('97^['^^2;,y,t])e‘ 


mt 


m=0 




E ^'^[b]l.fc(e^“)"'''" 'He|“’ ( 2 ^) 


Afterwards we note that [[l]].fc(e^'‘)""■’” " (see Appendix A, dSH)), 

■-[bU T -[bll 

E Ep (^1 E ^bAneri (x) 


J J =0 


m=0 




.7=0 


The sum by k represent the umbral composition (l57)l for the generalized Hermite polynomials HeJ“^ (x) [T^, which 
yields 

■-[bll 


m—0 


dx” 


j=o 
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As a result we obtain the following intermediate expression for the propagator 

^ 2M 2 'ijir ( \ 2 1 


i=i 


m—O 


We split the last expression in two terms = Ke + ^_r, where the first term contains error function in Eq{x^ y; t) 
whereas the second term contains elementary functions only 


-it 2M 

e 2 


ATb = ^(-l^e 4e y, t]He„[[j]] (x) 

^ ^ W{y) 


i=i 


— 2 W i i 2 1 rl-m 

Kb = E Zl, WK|-'.*.».'I)b=ob'’"'—(. r), 

m=l 




Changing operator L hy L with the aids of (l2^ we get 


-it 2M 

e 2 


ATs = —^ ^(-l)-^e 4 e j/, t]Heo-[[j]] (x) 


i=i 


VC(2/) 


-it 2M 

e 2 


Kb=—E(- lt 


£^Wjiy) ■ 


e 4 




W{y) 


■ 4 E 


(9j £i[J, X, y, t])j^g ^^^HeCT[[j]] (x). 


Let us consider the first term Ke- Though the function Eq contains non-elementary error function erf it can be seen 
that KE{x^y,t) is expressed in terms of elementary functions only. The error function will be cancelled due to the 
operator L as follows 




it 2M 


WM 


i=i 


-it 2M 

e 2 


E<-'> 




e 4 e 


i=i 


WM 

W{y) 


W{y) 

/ 

Eo[x,y,t]LJle^^^jYi i^) 

V 


2M \ 

E ^n,j{x)d^Eo{x,y\t) 

n—1 


W{x) 




where A„j (x) are some polynomials in x. By the definition (ED of the operator L we have LaE{e„\tjn (x) = 0, therefore 


-it y^-x^ 2M 


Ke = 


e 2 e 4 


2M-1 


2lT(y)lT(x) ^ 


^ A„^,{x)dr^Eoix,y-t) = 


n—0 


-it y^- 
0 2 6 4 


2M 


2M-1 


Koscix,y]t)e^ 4” e 2* ^(-1)^6 ^ A„,j(x) 

21T(i/)lT(x) ^ ^ 


n=0 


2i sinf 


^ Wi,„+i(x,y)e = 


i=o 


( i\ —it 2M— 1 y \ n n 

oscy, y; je --ia[b]]t^7 ) V- A„ ,(x) (- ) V u>/„+i(x, y)e “‘‘C 

21T(j/)TT(x) ^'v2isinf; E. 


=0 
—i t 


Thus we proved that Ke contains elementary functions only. Denote e by A. We can further simplify Ke writing 
it as a product of Kosc and a rational function of x, y and A, 


Ke = 


-^OSc(^; y■> ^') 

W{y)W{x) 


2A4 2Ad —1 / » \ 71 TL 

A^(-irA'^[W'ltj (?/) An,jix) ( ^ _ ^2 ) = 

j —1 n—0 ^ ^ l—O 
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4M+cr[[-l]]-l 
Kosc\^y 1 n—l 


W{y)W{x) 


(1- A2)2M-1 


where Qn{x,y) are some polynomials. 

Following the same line we will consider the second term. It can be seen that structures oi Ke and Kr coincide, 
and Kr reads 


Kr = 


4M+cr[[-l]]-l 

r. ( ..N E Qn{x,y)\^ 

-^osc(^5 y? n=0 


W{y)W{x) (l-A 2)2M-i 
Combining now Kr and Kr we obtain the following expression 


<t[[-1]] + 1+4M-2 

K (xvt) ^ 

J^(T _ J^OSC\X, y,l) „_Q 


^osciXf y, t) 


cr[[-l]] + l+4M-2 _ 

E Q^(a:,y)A^ 

n=0 


W{y)W{x) 


(1-A2) 


2M-1 


WT[hcrix),x]WT[ha{y),y] (1-A2) 


2M-1 


where some constant multiplier WT[hcr{x),x] = CW{x) is absorbed by redefinition of polynomials Q. 

Consider the limit t —?> 0, where K°' —>■ 5{x — y) and iFosc —t 5{x — y). From here it follows that the rational 
A-depending factor has no pole, therefore 


cr[[-l]] + l+4M-2 


E Qn{x,y)y' (i-A^) 


2M-1 


a[[-l]] + l 


(1-A2) 


2M-1 


E QZi^.v)^" 

,, _ -= E QU=:.«)y 

(1 A ) 


and 


<t[[-1]] + 1 

X! Qki^^y) = Wr[/i,^(a;),a;]Wr[/i<^(y),y]. 

k^O 

We finally get a rational anzatz for the propagators 

^OSC (^5 y ; 


K'^{x,y-t) = 


Wr [ha- (x ), a;] Wr [ha {y),y] 


a[[-i]]+i 

Qk{x,y)e~' 


kt 


(34) 


(35) 


fc =0 


where Q1{x,y) = Q1{y,x) are some polynomials to be determined. 


B. Rational anzatz for the propagators 

Substituting into the rational anzatz (1351) expansions of propagators in terms of eigen-functions 
^2,2°° 

Kosc{x,y-t) = e hn{x)hn{y)X^, A = e““, (36) 


n=0 


Ka{x,y,t) = \ ^ K{x)hl{y)X' 

m[ha{x),x\m[ha{y),y\ 


(37) 


we obtain a system of equations for the polynomial coefficients Q1{x,y). The solution of this system is given by a 
recursive procedure 


Qk{x,y) = iKix)K{y) - '^^Ql-jix,y)hj{x)hj{y)\ , 0 < A: < cr[[-l]]-f 1. 


(38) 













1. Nonlinear connection Lemma for the exeeptional Hermite polynomials 


Lemma. Given a Krein-Adler sequence a = + 1}) the corresponding family of (formally 

normalized) polynomials hf^{x) obeys the following relation 

<t[[-1]] + 1 

^ hm-k{x)hm-k(y)Ql{x,y) = hZ,{x)h^{y), (39) 

k=0 

where polynomials Q1 are given by (1551) . 

Proof. The proof follows from (1^ . (l36|) . (l37l) and (l38l) □ 

Relation (13911 can also be written in terms of wave functions 

<t[[-1]] + 1 

X) 'lfm-k{x)'ll}m-k{y)Q%{x,y) 

- = rmi^Wmiy) • ( 40 ) 

E Qlix,y) 

k=0 

There are following properties of the -polynomials: 

1) Symmetry 

Qkix,y) = Ql{y,x) 

2) Parity 

Qki-x,-y) = Ql{x,y) 

QllH-x,y) = (-l)^‘^®®''«^Q2fe(a:,y), Q 2 k+ii-x,y) = -{-l)^‘^‘'^^o'>Q 2 k+i{x,y) 

kl 

where deg/ij = E (o-[b']] - J + 1) “ W\- 
i=i 

In the Appendix B we present an example of a non-Krein-Adler sequence tr = {1} when the nonlinear connection 
lemma is also holds. We suppose that (1551) holds for an arbitrary a. 


2. x-Mehler formula 


We first recall the Mehler formula 


“ A” 1 ->.^(x^+v^)+2Xa:y 

VHe„(a:)He„(j/)— = _ e W^) 

Z—^ T7.' 

n—0 


n\ \/l - 


(41) 


Using the nonlinear connection lemma (1391) we obtain the following generalization of the Mehler formula to the case 
of exceptional Hermite polynomials 




-X^(x^+V^)+2Xxy 
e 2(i-A^) 


a[[-ll]+l 

Q‘"j(.x,y)\P 

i=o 


(42) 


3. Alternative form of the rationally extended Harmonic oscillators 
Consider integral kernels of Hamiltonian operators 

OO 

Hosc{x,y) = [-9E + Vosc{x)] S{x -y) = Y nifn{x)^n{y) , 


(43) 


n=0 


H'"{x,y) = [-dl^ + V^ix)] 5{x-y) = Y^^nix)'>Pniy) ■. 


n—0 


(44) 
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using dini) we can represent the second kernel as follows 


H'"{x,y) = 


a[[-l]] + l 

E kQ%{x,y) 


^xx “I" ^sc(a^) 


fc =0 


^[[-i]]+i 

E Q'j{x,y) 

j=o 


5{x-y). 


(45) 


From here is follows that 


^[[-i]]+i 

E kQl{x,x) 


AV^ix) = -2dxxlog[MVi[ipa{x),x]] = -■ 


k=0 


^[[-1]] + 1 

E Qj{x,x) 

j=0 


4- Green functions 

Consider the Green function (resolvent kernel of the Hamiltonian operator) 


CXJ 

G{x,y\E) = i j K{x,y;t)e'^*dt 


The transformation formula for the propagators implies also the following relation for the Green functions 

1 


Ga{x,y;E) = 


W[K{x), x]W[h^{y),y] 


a[[-l]] + l 

osc {x,y;E-k) 


(46) 


fc=0 


C. Examples 


Using first and second excited states of Harmonic oscillator 
ipi{x) = pixe tp 2 {x) = P 2 {x^ - l)e 

we obtain a perturbed Harmonic oscillator potential da 

V{i,2}N = ^ + 2(^1 + 2 ^^) 

Connection polynomials read 

p{i,2} _f 1 xy x'^y'^ + x'^ + y'^ - I xy) 
b/{o,i,2.3} - 2^’ Att ’27rJ' 

The propagator for the Schrodinger equation with Hamiltonian 
expression 




(47) 


-a2+u{i’2} 

(x) has the following compact 


(48) 


Next simple expression for the propagator we can obtain using second and third excited states of Harmonic oscillator 
tp 2 {x) = P 2 {x‘^ - l)e V'3(a;) = P3x{x‘^ - 3)e . 

In this case we obtain two-well perturbed Harmonic oscillator potential (potentials Vjfc fc+i} have k shallow minima 
at their bottom) 


^{2,3}[x] = Y + 2 


^1 -b 4x^ 


— 9 \ 
{x^ + 3)2 ) 


(49) 
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Connection polynomials read 

(a;^ + + 1) xy{3 — x^y'^) x^y‘^ + 3x^ + 3y^ — 12x'^y^ — 3 xy{x^y‘^ — 3) — 1) 


n{2.3} 

{0,1,2,3,4} ’ 0JJ. > 247r ’ Gtt ’ 47r 

The propagator for the Schrodinger equation with Hamiltonian _ _q 2 _|_ ^{2,3}^^;) reads 

^{ 2 , 3 }. (. 8ismt[xy{x^y^-3)-3{x^+y^)cost-3i{x^y^ + l)smt] \ 

(x,y,t)-e Kosc[x,y,t) (3 + ^ 4)(3 + y 4 ) j 


(50) 


III. CONCLUSIONS 


Propogators 

a[[-i]]+i 

E QU^,y)e-^^* 

K'^{x,y]t) = Kosc{x,y,t)—^^^--^^ - (51) 

E Qlix,y) 

fc=0 

present a new example of Feynman path integrals that can be calculated analytically [l^ . The key formula (l38l) 

Qk{x,y) = (^kix)K{y)-Y^^Ql_j{x,y)hj{x)hj{y^ , 0 < A: < cr[[-l]] + 1, (52) 

which define nonlinear connection between x-Hermite and Hermite polynomials can be easily realized in any computer 
algebra system. One can use these propagators to test various approximations for the corresponding path integrals. 
Using these propagators wave-packet dynamics in multi-well potentials can be studied analytically. 


Appendix A. Rescaled Hermite polynomials, Appell sequences and nmbral composition |l6l| 

The generating function of the probabilistic Hermite polynomials reads 
/ ^2\ °° 

R[z,x] = exp = X! Hen(x)— . 

n 

Heyj(x) = ^ ^ 7 

k^O 

Rescaled Hermite polynomials 

HeL“'(a;) = a^He™ , 

K'k = « " hn,k, 
with the following generating function 

R[zVa, -^] = exp (xz - ^ lle]^'^{x)^ . 

form an Appell sequence of polynomials. 

Let An(x] and Bn{x) be two Appel sequences of polynomials [l^ generated by functions S'[a;,g], R[x,g], 

oo „ 

s[x,g] = s(g)e“s = , 

n=0 


(53) 


(54) 

(55) 
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TL 

R[x,g] = r{g)e^o = ^ > 

n—0 


where An(x) = ^n,kX^ ,Bn{x) = J2 K,kX^■ 
fc =0 fe =0 


Define the umbral composition by the following formula 


n n k 

(-^n O -^)(^) ~ ^ ^ ^n,kBk{x) = EE 0'n,k^k,jX . 

k—0 k—0 j—0 


(56) 


Then the sequence Cn(x) = {An o B){x) also is an Appell sequence. Its generating function, denoted by U[x,g\ = 
S[x,g]o R[x,g\, has the following form 

°° n 

U\x,g] = s{g)r{g)e^3 = ^ C'„(x)^ . 


We can apply these facts to the sequences of generalized Hermite polynomials Hen\x) 


exp I xz 


az 


o exp I xz 




= exp I xz 


{a + I3)z'^ 


As a result we get the umbral composition of rescaled Hermite polynomials 

n 

o (x) = = He]^+^^{x). 




(57) 


Appendix B. Q-polynomials for non-Krein-Adler sequence a 


Here we consider an example of the polynomial connection when a = {1}, which is a non-Krein-Adler sequence 
[H, and hence, when polynomials hn do not represent a sequence of orthogonal polynomials. Nevertheless, by direct 
calculations we can verify that the Non-linear connection lemma is valid in this case. Consider formally normalized 
product of two wronskians 




1 

X hm{x) 

y hm{y) 

^/^{l — m) 

1 h'„,{x) 

1 h'nJ^y) 


(58) 


{xy/mhm-i{x) - hm{x)) [yy/mhm-iiv) - hmiy)) ■ 


^/^{l — m) 

Using the recurrence relation for the normalized probabilistic Hermite polynomials 


x\/mhm-i{x) = mhm{x) + m{m — l)hm- 2 {x) 


we get 




-\/^(l — m) 


\ 

-1 / 


hm{x) + 


\/{m-l) 


{^m - \]hm{x) + y m{m - l)hm- 2 {x)^ (^[m - l]hm{y) + \/m{m - l)hm -2 
hm- 2 {x^ {[m - l]hm{y) + \/m{m - l)hm -2 


[-hm{x)hm{y) + mhm{x)hm{y) + \/m{m- l){hm{x)hni-2{y) + hm{y)hm-2{x)) + mhjn-2{x)hm-2 

VZTT 

-^=hm{x)hm{y) - ^{xhm-1 (x) - y(TO - l)hn,-2(x))(yhni-i(y) - y(m- l)h„,_2(y)) 

yzir v2'7r 

y/{m-l) (^{xhm-l{x) - y(m - l)hm-2{x))hm-2{y) + iyhm^i{y) - \/{m - l)hm-2{y))hm-2ix)^ 
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1 


1 


mhm-2{x)hm-2{y) 


- hm{x)hm{y) - ■j={xhm-l{x)yhm-l{y) + (to - l)hm-2{x)hm-2{y)) 

V 27r y2iT 

+ — 2(m — m —2 (j/) 

V 27r 


1 


1 


mhm- 2 {x)h m—2 (y) 

{hm{x)hm{y) - Xyhm-l{x)hm-l{y) “ /lm-2 (a;)/lm-2 (?/) ) 


That is, polynomials (y) satisfy to the following 3-term representation 

2 

= '^hm-k(.x)hm-k{y)Q]^^{x,y) , 

where 


A;=0 




= - 


Q^^(.x,y) = - 


xy 

1 


We also verified by computer algebra that the nonlinear connection Lemma is valid for arbitrary cr with (t[[—1]] < 5. 
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